Development of a configuration-interaction + all-order method for atomic calculations 
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We develop a theoretical method within the framework of relativistic many-body theory to accu- 
rately treat correlation corrections in atoms with few valence electrons. This method combines the 
all-order approach currently used in precision calculations of properties of monovalent atoms with 
the configuration-interaction approach that is applicable for many-electron systems. The method 
is applied to Mg, Ca, Sr, Zn, Cd, Ba, and Hg to evaluate ionization energies and low-lying energy 
levels. 
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I. INTRODUCTION 

The development of the relativistic all-order method 
where all single and double excitations of the Dirac- 
Hartree-Fock wave function are included to all orders of 
perturbation theory led to accurate predictions for ener- 
gies, transition amplitudes, hyperfine constants, polariz- 
abilities, C3 and Ce coefficients, isotope shifts, and other 
properties of monovalent atoms, as well as the calcula- 
tion of parity-violating amplitudes in Cs, Fr, and Ra^ 
[l|, [3, H, 0, S @1 ■ This method was also used to calcu- 
late magic wavelengths Q as well as black-body radiation 
shifts [8| and quadrupole moments [9] that are of inter- 
est to atomic-clock research. The all-order method is 
designed to treat core-core and core- valence correlations 
with high accuracy. It is one of the most accurate meth- 
ods currently being used in the atomic structure calcula- 
tion. However, its applications so far have been limited 
to monovalent systems. Readers are referred to Ref. 
and references therein for a review of this method and its 
applications. 

Precision calculations for atoms with several valence 
electrons require an accurate treatment of the very 
strong valence-valence correlation; a perturbative ap- 
proach leads to significant difficulties. The complexity of 
the all-order formalism for matrix elements also increases 
drastically as the number of valence electrons increases; 
for example, the expression for all-order matrix elements 
in divalent systems contains several hundred terms in- 
stead of the twenty terms in the corresponding mono- 



*URL: |http: //wMw.udel . edu/~msaf rono[ 
lmsafrono@udel.edul 



Electronic address: 



valent expression. Therefore, we found it impractical to 
develop a direct extension of the all-order approach to 
more complex systems, both due to the large valence- 
valence correlation corrections and the very large number 
of terms noted above in the matrix element formulas. 

A more promising method for the study of atomic 
properties of more complicated systems that combined 
configuration interaction (CI) and many-body perturba- 
tion theory (MBPT) was developed in Ref. [Si- The 
CI-I-MBPT method was appfied to the calculation of 
atoniic properties of various systems in a number of works 
(see [ul [il [ll H [11, [ii, [13, [11, [li, [13] and references 
therein) and to the calculation of PNC amplitudes in 
Tl and Yb [H], [H]. The strengths of the configuration- 
interaction method are broad applicability and all-order 
treatment of the valence-valence correlation corrections. 
However, the precision of the CI method is generally dras- 
tically limited for large systems by the number of the 
configurations that can be included. As a result, core 
excitations are neglected or only a small number of them 
are included, leading to a significant loss of accuracy for 
heavier systems. The CI-I-MBPT approach allows one to 
incorporate core excitations in the CI method by con- 
structing an effective Hamiltonian that incorporates 
certain perturbation theory terms. The CI method is 
then applied to the modified H'^^ to obtain improved en- 
ergies and wave functions. 

Because of the rapid increase in the number of terms 
of the MBPT expansion, the CI + MBPT approach be- 
comes impractical already in the third order of MBPT. 
For that reason, the CI + MBPT approach is usually re- 
stricted to second order. Some higher-order corrections 
can be accounted for by introducing screening coefficients 
to second order diagrams. These screening coefficients 
can be found either by averaging two-electron second or- 
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der diagrams, or by semi-empirical fitting of experimental II. 
energies. The second-order expression for one-body cor- 
rection to the Hamiltonian is corrected by the all-order 
chains of such terms in some works (see, for example, 
Ref. [231). In 2004, a modification of the effective Hamil- 
tonian using the all-order pair equations was proposed 
and tested on a "toy" 4-electron model [24*] . An efficient 
method of including core-valence correlations into the 
configuration interaction (CI) calculations was presented 
by Dzuba and Flambaum in [2^. The CI Hamiltonian 
for N valence electrons was calculated using orbitals in 
the complete potential (the mean field produced by 
all electrons); the one- and two-body corrections to the 
effective Hamiltonian were obtained by using many-body 
perturbation theory with dominating classes of diagrams 
included to all orders. 

In the present work, we combine the all-order method, 
currently used in precision calculations of properties of 
monovalent atom, with the configuration interaction (CI) 
approach. In the CI -I- all-order approach, core excita- 
tions are incorporated in the CI method by constructing 
an effective Hamiltonian using fully converged all-order 
excitation coefficients. Therefore, the core-core and core- 
valence sectors of the correlation corrections for systems 
with few valence electrons will be treated with the same 
accuracy as in the all-order approach for monovalent 
atoms. Then, the CI method is used to treat valence- 
valence correlations. This method is expected to yield 
accurate wave functions for subsequent calculations of 
atomic properties such as lifetimes, polarizabilities, hy- 
perfine constants, parity-violating amplitudes, etc. The 
present work is motivated by the urgent need for pre- 
cision calculations of atomic properties of heavy atoms 
with few valence electrons for applications such as atomic 
clock research, quantum information, study of fundamen- 
tal symmetries, searches for variation of the fundamental 
constants, and tests of high-precision experimental meth- 
ods. The development of the CI -t- all-order method is 
also aimed at filling the long-standing gap between the 
accuracy of theoretical and experimental parity- violation 
studies in systems with few valence electrons. Atomic 
properties of various atoms and ions are also of interest 
for astrophysics applications. 

Our method is generally applicable, i.e. not restricted 
to the specific type of the system. We test the present 
approach on the calculation of the energy levels of Mg, 
Ca, Sr, Cd, Zn, Ba, and Hg to demonstrate a signifi- 
cant improvement in comparison with CI -I- second-order 
MBPT values. We also discuss calculations of transition 
matrix elements and polarizabilities. 

We provide a brief description of the all-order and CI l*^*^) 
-I- MBPT formalisms in Sections HIl and Hill respectively. 
The CI + all-order approach is described in Section [TVl 
Finally, we present results for removal energies in diva- 
lent systems calculated using the CI -|- all-order approach 
and discuss perspectives for further applications and de- 
velopment. 



RELATIVISTIC ALL-ORDER METHOD FOR 
MONOVALENT SYSTEMS 



Our point of departure is the relativistic no-pair Hamil- 
tonian H = Hq + Vi obtained from QED by Brown and 
Ravenhall [Hi: 
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Here, Vijki are two-particle matrix elements of the 
Coulomb interaction gijki, or Coulomb + Breit interac- 



tion gijki+hjku and Vhf = Y.a 



v^aaj) \s frozen- 



core Dirac-Fock potential. The summation index a in 
Vhf ranges over states in the closed core. The quan- 
tity Ei in Eq. ([1]) is the eigenvalue of the Dirac equation 
h{r)(pi{r) = ei4>i{r), where 



h{r) — COL ■ p + [3m(? 



Z 
r 



U{r) 



(3) 



In our previous all-order calculations of monovalent 
atoms, we took U to be frozen-core y^-i potential, 
U — Vhf- Such a choice greatly simplifies the calcu- 
lations since the second term in Eq. ^ vanishes in this 
case. In this work, we use the same type of potential 
(]/^~2 for divalent systems) , but different potentials 
may be used in the future. 

In the coupled-cluster method, the exact many-body 
wave function is represented in the form (27| 
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exp(5)|vI/(")), 



(4) 



where is the lowest-order atomic state vector. The 

operator S for an N-electron atom consists of "clus- 
ter" contributions from one-electron, two-electron, • ■ • , 
N-electron excitations of the lowest-order state vector 
|*(0)): S^Si + S2 + --- + Sn. 

The all-order method described in detail in Rcfs. 
is a linearized version of the coupled-cluster method, 
where all non-linear terms in the expansion of the ex- 
ponential are omitted; the all-order wave function takes 
the form 



{1 + S-l + ^2 + ^3 



5Ar}|*(0)). (5) 



Restricting the sum in Eq. ([5]) to single, double, and 
valence triple excitations yields the following expansion 
for the state vector of a monovalent atom in state v: 



(6) 
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where the indices m, n, and r range over all possible 
virtual states while indices a and b range over all occupied 
core states. The lowest-order wave function l^-i"^} is 



(7) 



where |5'c) is the lowest-order frozen-core wave function. 
The quantities Pma, Pmv are single-excitation coefficients 
for core and valence electrons; Pmnab and Pmnva are core 
and valence double-excitation coefficients, respectively; 
Pmnrvab are the valence triple excitation coefficients. In 
the single-double (SD) implementation of the all-order 
method, only single and double excitations are included. 
In the (single, double, partial triple) SDpT variant of the 
all-order method, valence triple excitations are included 
perturbatively as described in Ref. 0]. 

To derive equations for the excitation coefhcients, the 
state vector IvP^) is substituted into the many-body 
Schrodinger equation iJj^'t,} = £'|^t,), and terms on 
the left- and right-hand sides are matched, based on 
the number and type of operators they contain, leading 
to the equations for the excitation coefficients given in 
Safronova et al. • 

We note that all non-linear terms at the single-double 
level have been added in the formulation of the all-order 
method in Ref. [2^ . This version of the all-order method 
is equivalent to the coupled-cluster single-double (CCSD) 
method for a finite basis set. It was shown in Refs. (29l.[30| 
that both non-linear terms and complete valence triple 
excitations must be included to improve the accuracy 
of the linearized coupled-cluster SD method. In the 
present work, we use linearized SD variant of the all-order 
method since the SD method already leads to excellent 
results for a large number of the atomic properties and 
is computationally efficient. 

The resulting SD all-order equations for valence exci- 
tation coefficients are: 



^ ^ QmhnrPnrvb 
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^ ^ QmnrsPrsvb 
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QmnrbPrv 



9cnvbPmc 



^ ^ QcnrbPmrvc 
rc 



V ^ b 
m ^ n 



where SE^ is the valence correlation energy SE-^j — Ey — 
; ^ij — H" ^ J ; and Pmnvb — Pmnvb Pnrnvb- The Corre- 
lation correction to the energy of the state i;, is given in 
terms of the excitation coefficients by 



6Ey — ^ gya 



2Pma ~(~ ^ ^ QabvmPmvab 
mab 



^ ^ QvbmnPmnvb- 
mna 

(10) 



Equations for core excitation coefRcients pma and Pmnab 
are obtained from the above equations by replacing the 
valence index by a core index a and removing SEy from 
the left-hand side of the equations. We note that the 
right-hand side of the valence energy equation is identical 
to the right-hand side of the equation for pmv with m = v. 
Equations for the correlation energy and all excitation 
coefficients are solved iteratively. Every iteration picks 
up correlation terms corresponding to the next higher 
order of perturbation theory until the correlation energy 
converges to sufficient numerical accuracy. Therefore, the 
all-order approach includes dominant MBPT terms to 
all-orders. 

Matrix elements for any one-body operator Z = 
X^i" ctiO'i are obtained within the framework of the 



all-order method as 



(11) 



where j^I'i,) and |^tu) are given by the expansion ([6]). In 
the SD approximation, the resulting expression for the 
numerator of Eq. consists of the sum of the DF ma- 
trix element Zwv and 20 other terms that are linear or 
quadratic functions of the excitation coefficients. The ad- 
vantage of this approach is that the expression in Eq. (jlip 
does not depend on the nature of the operator Z, only 
on its rank and parity. Therefore, matrix elements of 
any one-body operator may be calculated with the same 
general code. 

The complexity of the all-order formalism for matrix 
elements increases drastically with the number of valence 
electrons. We have derived the expression for all-order 
matrix elements in divalent systems; it contains several 
hundred terms instead of the twenty terms in the cor- 
responding monovalent expression. Therefore, it is im- 
practical to extend the all-order method to the case of 
more complicated atoms in its present implementation 
directly, i.e. to start single-double expansion from the di- 
valent lowest-order wave function containing a^^all'^c)- 
Such an approach also leads to "intruder state" problems 
well-known in the perturbation expansions based on the 
Rayleigh-Schrodinger implementation of the MBPT. 

We note that the relativistic couple-cluster method has 
been successfully applied to calculation of the energies 
and electron affinities in systems with few electrons (see 
Refs.[3l|, m, m, m, ll^ and references therein). It is 
rather difficult to apply this method for calculation of 
other atomic properties, such as transition matrix ele- 
ments. However, using the finite field technique one can 
use this method, for example, to calculate quadrupole 
hyperfine constants in such a heavy atoms as Au [36|. 



III. CI + MBPT METHOD 

A combination of the configuration-interaction (CI) 
method and perturbation theory was developed in 
Ref. It was based on the Dirac-Fock code [s^l and 
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CI code [ssj . This approach has been later applied to the 
calculation of atomic pro perties of various systems in a 
number of works (see |ll|, M, M,M^ M,M, d and 
references therein). In Refs. HH, HI the CI + MBPT 
method was used to calculate PNC amplitudes in Tl 
and Yb respectively (in the latter case the nuclear-spin- 
dependent amplitude was calculated). 

In the CI method, the many-electron wave function is 
obtained as a linear combination of all distinct states of 
a given angular momentum J and parity 



(12) 



in other words, a linear combination of Slater determi- 
nants of proper symmetry from a model subspace [lO] . 

Energies and wave functions of low-lying states are de- 
termined by diagonalizing an effective Hamiltonian: 



(13) 



where Hi represents the one-body part of the Hamil- 
tonian, and H2 represents the two-body part (Coulomb 
or Coulomb + Breit matrix elements Vijki)- We use 
Coulomb matrix elements gtjki in the present work. The 
resulting wave functions are used to calculate matrix 
elements and other properties such as polarizabilities, 
parity-violating amplitude, etc. The precision of the 
configuration-interaction method is drastically limited 
for large systems by the number of the configurations 
that can be included. Consequently, core excitations are 
entirely omitted or only a small number are included, 
leading to a significant loss of accuracy for heavy atoms. 

The CI 4- MBPT approach allows one to incorporate 
core excitations in the CI method by including certain 
higher-order terms in an effective Hamiltonian (jl3p . The 
CI method is then applied as usual with the modified H'^^ 
to obtain improved energies and wave functions. Some- 
what different versions of the CI-I-MBPT method exist; 
here, we describe the approach used in the present work 
and follow the designations of Ref. . 

In the CH-MBPT approach, the one-body part Hi is 
modified to include the correlation potential Ei that ac- 
counts for part of the core- valence correlations. 



Hi — > Hi + Si. 



(14) 



(2) 

Either the second-order expression, : o^' all-order 
chains of such terms can be used (see, for example, 
Ref. [2^). The latter approach corresponds to replacing 

Dirac-Fock orbitals by Brueckner orbitals. The second- 

(2) 

order matrix elements (E^ )yx are given by 



Qmyab Qmxab 



niao 



rnnxa 9mnya 



(15) 

We use the same designations as in Section [Hi indices 
from the middle of the alphabet m and n range over all 
possible virtual states while indices a and b range over 



all occupied core states. The one-particle energies Ci are 
written together as — ei + Cj for brevity. The sum- 
mation over index i implies the sum over the quantum 
numbers rii ki rrii . 

The CI+MBPT approach is based on the Brilloiun- 
Wigner variant of MBPT, rather than the Rayleigh- 
Schrodinger variant. Use of the Rayleigh-Schrodinger 
MBPT for systems with more than one valence elec- 
tron leads to a non-symmetric effective Hamiltonian and 
to the problem of "intruder states". In the Brilloiun- 
Wigner variant, the effective Hamiltonian is symmetric 
and accidentally small denominators do not arise; how- 
ever. El and S2 became energy dependent. Specifically, 
the one-body correction Ei depends on the energy iy (see 
Eq. Ideally, the energy iy should be calculated 

from the particular eigenvalue of the effective Hamilto- 
nian (|13p . In practice, we use several approaches. The 
simplest and the most practical one is to set the energy 
iy to the Dirac-Fock energy of the lowest orbital for the 
particular partial wave. For example, we use ins = ess for 
all ns orbitals of Mg system. This approximation usually 
works reasonably well for atomic states belonging to the 
lowest configurations of a given symmetry. Another ap- 
proach is to set the energy of all orbitals for a particular 
partial wave to a certain value, one value is specified for 
each partial wave. This approach allows one to generate 
better atomic wave functions for subsequent evaluation of 
the atomic properties by selecting the values of e so final 
energy eigenstates are tuned to the experimental values. 
We have also developed a more elaborate method that 
involves calculating derivatives of the Ei and E2 with re- 
spect to e that allows to adjust the effective Hamiltonian 
as suggested in (T^]. Our implementation of the Cl+all- 
order method permits us to use any of these strategies. 

Performing analytical sums over all magnetic quantum 
numbers yields the expression 



(Ef)) 



(2)\ _ 



yx 



EE 

mab K 

EE 

mna K 



1 Xj^ [ray ah) (mxab) 



m [jy 



1 XK{fnnxa)ZK{'mnya) 



[K] [jy 



, (16) 



where K is multipolarity restricted by conventional tri- 
angular rules, [K] = 2K + 1, and the summations over 
all lower-case indexes i now designate sums over nt and 
Ki. We use similar designations for the sums listed below 
to avoid explicitly writing out all quantum numbers in 
all sums. The quantity Xximnab) is 



Xnimnab) = 



(17) 



(-1)"^ (Km ||C^ II Ka) («„ ||C^|| ^b) RK{mnab), 

where RKimnab) are (relativistic) Slater integrals and 
(kto ||C^|| I'^a) are reduced matrix elements of a normal- 
ized spherical harmonics. ZKimnab) is given by 



ZxiTnuab) = XniiTinab) 



(18) 



EM 

K' 
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The two-body Coulomb interaction term H2 is modi- 
fied by including the two-body part of the core-valence 
interaction that represents screening of the Coulomb in- 
teraction by valence electrons; 



(19) 



where E2 is calculated in second-order MBPT in 
CI-I-MBPT approach. The second-order matrix elements 
are given by 



Qvwcd Qmncd 



(20) 



Qwrnc Qmrvc 



m 4^ n 
V w 



Performing an angular reduction leads to 



(21) 



jm jv ^ l J ^ 

L K' \ L K' jd 



cd LK' 

Xl (vwcd) Xk' (cdmn) 



(^^ly^u+jn+K ZK{wrnc) ZK{iTLrvc) 



E 

rc 

E 



(^_lp\;+jm+A ZK(vrmc) ZK{nrwc) 



[K] 



MBPT corrections associated with terms Ei in Eq. p4)) 
and E2 in Eq. (fT9|) typically grow with nuclear charge 
Z leading to a deterioration of the accuracy of the CI 
-I- second-order MBPT results for heavier, more com- 
plicated systems. The order-by-order extension of this 
method does not look promising for two reasons. First, 
the complexity of the MBPT expansion for systems with 
more than one valence electron already makes third-order 
calculations impractical. Second, the convergence of the 
MBPT series is not well studied, but it is known that 
third order is often less accurate than second order. This 
is why it was so important to develop an all-order exten- 
sion of the MBPT method for monovalent systems. 



IV. CI + ALL-ORDER METHOD 

In the CI + all-order approach, corrections to the ef- 
fective Hamiltonian Ei and E2 are calculated using the 
all-order method, in which the effective Hamiltonian con- 
tains dominant core and core-valence correlation correc- 
tions to all orders, as discussed in Section |TT1 The core- 
core and core-valence sectors of the correlation correc- 
tions for systems with few valence electrons are treated 
in the all-order method with the same accuracy as in the 
all-order approach for the monovalent systems. The CI 
method is then used to evaluate valence-valence correla- 
tions. 



First, we express the all-order equations Eqs. (|81 fTC)l) 
in terms of matrix elements of Ei and E2 and explicitly 
include the energy dependence. We also need to add an 
all-order equation for the excitation coefficients Pmnvw to 
obtain E2. This equation is equivalent to Eq. ([9]) with 
core index b replaced by valence index w. 

El and E2 are essentially the all-order excitation coef- 
ficients Pmv and Pmnvw ■ 



^ma — Pma {^a ^m) 

'-'mnab — Pmnah {^ab ^mn) 

-'mnva — Pmnva {^v ^~ ^mn) 

Emtj — (El)^^^ Pmv {^v dm) 

^mnvw — {"^2) mnvw Pmnvw {dy -\- ^mn) 



The quantities E„ia, I^mnab, ^mnva are used in the all- 
order iteration procedure but do not explicitly appear in 
the effective Hamiltonian. The core equations for pma 
and Pmnab are not modified from the original all-order 
monovalent code. The excitation coefficients pma and 
Pmnab are simply multiplied by the appropriate energy 
differences to obtain the terms E^a and ^mnab needed by 
other programs. Re- writing the other all-order equations 
in terms of E and removing terms that will be otherwise 
double-counted by the CI part of the calculations, we 
obtain the following set of equations: 



(Si) 
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7ib 



9inbv7i ^7ih 



(22) 
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bcTi 



97nbnr ^nrvb 
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^ ^ 9cdvb ^mncd 

, ^cd ^mn ~h C-y €y 
cd 



(23) 
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QTTlCvb Ej7,C 
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9cnrb Er 
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Ci; + ec - er^ 



E 



^Tncvr ^ 



rncb 



^cb ^rir ~l~ 



E97ncvr ^nrcb QcTivr ^7nrcb 

_ zr„„ _|_ _ z:,. Z-^ 



^cb ^Tir ~l~ ^cb ^mr 
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(24) 



Qcdvw ^rancd 



Qcnvw ^mc 

















Qmcvw ^Tic 




fc 








Qcnrw ^mrvc 




El. 








Qcnrw ^rrnvc 




eu 








Qincrw ^mvc 




El. 




~ 




9mcvr ^mcw 








- 




Qmcvr '^rnwc 




ec 




- ev 




Qcnvr ^rrawc 








- £v 



E 

cd 

E 

c 

E 

c 

E 

cr 

E 

cr 

E 

cr 

E 

cr 

E 

cr 

E 



The energy denominators are now explicitly written out 
and the energy dependence is introduced following the 
prescription of the CI + second-order MBPT approach. 
Putting e„ = e^, yields the original all-order equations ([5]) 
and ^ for monovalent systems up to SEy terms on the 
left-hand side. The term containing prv is removed from 
Eq. since it is included in the CI calculation. Eq. (IM)) 
for 'Emnvw does not have terms that would corresponds 
to the first, third, and fifth terms of Eq. for J^mnva, 
since these terms are accounted for by the CI as well. 
Therefore, no iteration is required in the last equation 
since T,mnvw does not appear on its right-hand side; the 
last equation needs to be evaluated only once after all 
other equations have converged. The last equation is also 
significantly faster to evaluate than similar equations for 
Smnab and Smnua sincc it docs not contain the term with 
four indexes over the virtual orbitals (term 3 in Eq. I23p . 
We note that this is the only equation that is not present 
in any form in the all-order code for monovalent systems. 

Below, we outline a step-by-step application of the CI 
+ all-order method. 

Step 1: A finite basis set is generated in a spherical 
cavity of radius R. All calculations are carried out with 
same basis set. The second-order MBPT is used to gen- 
erate (S^^^)t,u; and {'S'^'')mnvw, where m,n,v,w are any 
basis set functions. Highly-excited orbitals can be omit- 
ted without loss of accuracy. A subset of the basis set 
orbitals is selected for which all-order values of 
and {T.2)mnvw are to be obtained. In our calculation this 
set includes the first three ns, 'T-P3/2: ^'^3/2j '^'^5/2 

states for each system. We found that inclusion of the ad- 



ditional orbitals did not significantly change the results 
and treating the remaining corrections to the effective 
Hamiltonian in second order should be sufficiently accu- 
rate. This issue is addressed in more detail in the next 
section. 

Step 2: The all-order core pma and Pmnab excitation 
coefficients are obtained by the iterative solution of the 
corresponding equations in the appropriate potential (for 
example, V^~'^ for divalent systems) in the same way as 
for monovalent systems. The core correlation energy is 
used as a convergence parameter and is generally required 
to converge to 10~^ relative accuracy. The core excita- 
tion coefficients are multiplied by the appropriate denom- 
inators as described above to obtain Tima and Emnafc after 
the iterations are complete. 



Step 3: The core quantities Ema and ^mnab are used 
to obtain J^mv and J^mnva, again by an iteration proce- 
dure, for a large number of excited m,n, and v orbitals. 
The valence correlation energy for the state v is used as 
a convergence parameter. These steps are carried out in 
exactly the same way as our present monovalent all-order 
calculations with the omission of the valence-valence di- 
agrams as described above. The iterations of excitation 
coefficients result in the summation of the relevant classes 
of MBPT terms to all orders. We note that the term 
gives the all-order correction to the one-body part of the 
effective Hamiltonian. 

Step 4: The all-order expression for {'S2)mnvw cor- 
rections to the effective Hamiltonian are calculated using 
Eq. (|24p with previously stored, fully converged, values 

of ^ma; ^jnnah^ and S^^^^^. 

Step 5: CI calculations are carried out to generate 
accurate wave functions with the effective Hamiltonian 
constructed using Ei and S2 obtained in the previous 
steps. 

Step 6: The resulting wave functions are used to ob- 
tain various matrix elements and derived quantities such 
as PNC amplitudes. In the current CI + MBPT ap- 
proach [l6| , matrix elements are calculated by replacing 
"bare" matrix elements by the "dressed" matrix elements 
using the random-phase approximations (RPA). In this 
work, we use the same approach. This issue will be fur- 
ther discussed in the later section. 

The method described above treats electronic correla- 
tion in systems with several valence electrons in a sig- 
nificantly more complete way than the CI -I- MBPT ap- 
proach owing to the inclusion of the additional classes 
of MBPT terms in Si and addition of all-order (rather 
than second-order) correction in £2- We note that our 
present all-order code is capable of efficiently evaluating 
the large number of the core-valence all-order excitation 
coefficients needed for the implementation of the CI -I- 
all-order approach. 
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TABLE L Comparison of the CI, CI+MBPT and Cl+all-order ab initio results for the two-electron binding energies of Mg, 
Ca, Zn, Sr, Cd, Ba, and Hg with experiment. The energies are given in cm"'. The relative difference with experimental values 
is given in the last three columns in %. 

Element State Energies Differences with experiment (%) 







Expt. 


CI 


CI+MBPT 


Cl+all-order 


CI 


CI+MBPT 


Cl+all-order 


Mg 




182939 


179537 


182717 


182877 


1.86 


0.12 


0.03 


Ca 


4s' 'So 


145058 


139068 


145985 


145517 


4.13 


-0.64 


-0.32 


Zn 


4s" So 


220662 


204083 


218521 


219442 


7.51 


0.97 


0.55 


Sr 


5s' 'So 


134896 


127858 


136082 


135322 


5.22 


-0.88 


-0.32 


Cd 


5s' 'So 


208915 


188884 


210716 


208620 


9.59 


-0.86 


0.14 


Ba 


6s' 'So 


122721 


114898 


124956 


123363 


6.37 


-1.82 


-0.52 


Hg 


6s' ' So 


235469 


207652 


241152 


236626 


11.81 


-2.41 


-0.49 



V. RESULTS AND DISCUSSIONS 

We compare the results of our CI, CI + MBPT, and 
CI + all-order ah initio calculations for the two-electron 
binding energies of Mg, Ca, Zn, Sr, Cd, Ba, and Hg with 
experiment in Table [H Results for the energies of Mg, 
Ca, Cd, and Ba, counted from the ground state, are 
compared with experiment in Table [Til The same des- 
ignations are used in both tables. The energy values are 
given in cm~^. Relative differences of our results with 
experiment are given in the last three columns of Tables 
HI and Hn to illustrate the accuracy of each approach. 

The same parameters and basis set are used in all three 
calculations for each system. The finite basis set of 245 
orbitals that include I = ... 5 partial waves is formed 
in the spherical cavity with a 50 a.u. radius. The CI 
calculation includes only valence shell excitations as de- 
scribed above, i.e. CI calculation is carried out in the 
same way for all three cases. All summations over the 
excited states in the second-order and the all-order cal- 
culations are always carried out over the entire basis set. 

A sufficiently large number of the effective Hamilto- 
nian matrix elements are modified in the CI + second- 
order MBPT calculation. There is no need to include 
corrections to the entire Hamiltonian as the corrections 
from the remaining terms are negligible as described be- 
low. The Cl+all-order calculations include replacement 
of the most important Si and S2 terms by their all-order 
values. The remaining corrections from the effective 
Hamiltonian retain their second-order values as described 
above. We find that it is sufhcient to carry out all-order 
calculations for the first three ns, np\i2-, '^■^3/21 

states and modify the corresponding Ei and £2- 
For example, indexes m,n,v,w in (Si)mi> and {Ti2)mnvw 
in Ca calculation include 4s, 5s, 6s, 5pi/2, 6pi/2, 

4p3/2, 5p3/2, 6p3/2, 3(^3/2, 4(^3/2, 5^3/2, 8^5/2, 4^5/2, and 
5ci5/2 states. To test that the above number of the cor- 
rected Hamiltonian matrix elements is sufficient, we have 
carried out the following test in Ca: 

(1) the number of second-order Ei and E2 matrix ele- 
ments included was increased from 535 to 878 and from 
4 879 832 to 19 236 743, respectively; 

(2) the number of second-order Ei and E2 matrix ele- 



ments replaced by all-order values was increased from 168 
to 305 and from 592 634 to 2 898 122, respectively, by 
including the 7s, 7pi/2, Ip^/i, 6^3/2, 6^5/2, 4/5/2, 4/7/2, 
5/5/2, 5/7/2 states into the m,n,v,w index set. 

The ionization potential and most of the Ca levels 
that we considered shifted by less than 1 cm~^. The 
largest changes, observed for the 4s3(i levels, were still 
very small, 0.1 %. This was expected since the nd levels 
are known to be the most affected by the partial wave 
restrictions. We note that the energies and in 
Eqs. (|22ll24p were set in the present calculation to the 
corresponding Dirac-Fock values for the lowest state for 
each partial wave. The second-order calculations were 
carried out in the same way. Therefore, the results listed 
in Tables HI and [Til are completely ab initio. 

We find that the all-order ionization potential results 
are in significantly better agreement with experiment in 
comparison with the CI+MBPT values even in the case 
of Mg where the agreement with experiment is already 
excellent in the CI+MBPT approach. We also find al- 
most no deterioration in the accuracy of the two-electron 
binding energies from Ca to Hg; the all-order method 
reduces the differences with experiment by about factor 
of three in comparison with the second-order data. Sim- 
ilar improvements are observed for most of the excited 
states listed in Table HIl with the exception of 5fi6s states 
of Ba. The accuracy of the SD all-order approach is ex- 
pected to be lower for Ba since the SD method omits 
certain parts of the third-order energy correction associ- 
ated with valence triple excitations (last term in Eq. ([6])). 
This contribution was found to increase significantly for 
heavier alkalis [2]. The problem is corrected in the SD 
all-order method by explicitly adding the missing part 
of the third-order correction. Within our approach, this 
issue may be treated in an ah initio way by adding the 
valence triple excitation terms perturbatively to the all- 
order as was done for the monovalent systems in Ref. Q 
and removing terms that are accounted for by the CI. 
We note that factor of three improvement in the relative 
differences with experiment is still observed for the 6s6p 
energies calculated with the all-order method. 

We observe that essentially all of the states in Ba listed 
in Table nil are shifted by the same relative amount in the 
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TABLE II: Comparison of the CI, CI+MBPT and Cl+all-order ab initio results for the energy levels of Mg, Ca, Cd, and Ba 
with experiment. Two-electron binding energies are given in the first row for each element, the other values are counted from 
the ground state energy. The energies are given in cm~^. The relative difference with experimental values is given in the last 
three columns in %. 



Element State Energies Differences with experiment (%) 





Expt. 


CI 


CI+MBPT 


Cl+all-order 


CI 


CI+MBPT 


Cl+all-order 


3s^ ''So 


T O O 

182939 


179537 


182717 


182877 


1.86 


0.12 


0.03 




4II97 


40409 


41132 


41175 


1.91 


0.16 


0.05 


3s4s ^Sn 


43503 


42689 


43459 


43502 


1.87 


0.10 


0.00 




46403 


45119 


46318 


46384 


2.77 


0.18 


0.04 


3s3d^Di 


47957 


46972 


47892 


47936 


2.05 


0.14 


0.04 




47957 


46972 


47892 


47936 


2.05 


0.14 


0.04 


3s3(7 3 Do 


47957 


46972 


47892 


47936 


2.05 


0.14 


0.04 


:/ • > Lj J. y 


21850 


20906 


21780 


21833 


4.32 


0.32 


0.08 


3s3p^Pi 


21870 


20926 


21801 


21852 


4.32 


0.32 


0.08 


3s3p^P2 


21911 


20967 


21844 


21897 


4.31 


0.30 


0.07 


3s3p^Pi 


35051 


34488 


35053 


35065 


1.61 


0.00 


-0.04 




145058 


139068 


145985 


145517 


4.13 


-0.64 


-0.32 


3d4s^Di 


20335 


24200 


19927 


20335 


-19.00 


2.01 


0.00 


ZdAs^Do 


20349 


24201 


19949 


20355 


-18.93 


1.97 


-0.03 


3d4s^D3 


20371 


24203 


19982 


20386 


-18.81 


1.91 


-0.07 


3d4s ^ Do 


21850 


23853 


21620 


21965 


-9.17 


1.05 


-0.53 




31539 


30147 


31765 


31694 


4.42 


-0.72 


-0.49 


4s5s ^Sn 


33317 


31893 


33552 


33466 


4.27 


-0.70 


-0.45 




15158 


13509 


15474 


15338 


10.88 


-2.08 


-1.19 


4s4p^Pi 


15210 


13557 


15528 


15385 


10.87 


-2.09 


-1.15 


4s4p^P2 


15316 


13655 


15638 


15498 


10.85 


-2.10 


-1.19 


4s4p ^Pi 


23652 


23052 


23771 


23729 


2.54 


-0.50 


-0.32 


tr 2 1 o 
5s So 


208915 


188884 


210716 


208620 


9.59 


-0.86 


0.14 




51484 


44027 


51916 


51395 


14.48 


-0.84 


0.17 


5s6s "^5*0 


53310 


46153 


53788 


53272 


13.43 


-0.90 


0.07 




59220 


50634 


59697 


59015 


14.50 


-0.81 


0.35 


5s5d^_Di 


59486 


51292 


59881 


59259 


13.77 


-0.66 


0.38 




59498 


51303 


59893 


59271 


13.77 


-0.66 


0.38 




59516 


51320 


59911 


59291 


13.77 


-0.66 


0.38 


5s5p^Po 


30114 


24417 


30903 


30141 


18.92 


-2.62 


-0.09 


5s5p^Pi 


30656 


24875 


31451 


30646 


18.86 


-2.59 


0.03 


5s5p^P2 


31827 


25833 


32656 


31838 


18.83 


-2.60 


-0.03 


5s5p^Pi 


43692 


38902 


43970 


43607 


10.96 


-0.64 


0.20 




122721 


114898 


124956 


123363 


6.37 


-1.82 


-0.52 




9034 


11524 


9276 


9249 


-27.57 


-2.67 


-2.38 




9216 


11603 


9489 


9441 


-25.91 


-2.97 


-2.45 




9597 


11780 


9941 


9840 


-22.75 


-3.59 


-2.54 


6s5d^D2 


11395 


12753 


11878 


11727 


-11.92 


-4.24 


-2.91 


6s6p^Po 


12266 


9938 


13112 


12556 


18.98 


-6.90 


-2.36 


6s6p^Pi 


12637 


10269 


13484 


12919 


18.73 


-6.70 


-2.23 


6s6p^P2 


13515 


11010 


14391 


13819 


18.53 


-6.48 


-2.25 


6s6p^Pi 


18060 


16908 


18621 


18292 


6.38 


-3.11 


-1.28 



all-order approach, unlike the case of CI + MBPT. In this 
case, it is possible to carry out another calculation with 
different and resulting in final energies in very close 
agreement with experiment in order to get an improved 
representation of the wave functions for subsequent eval- 
uation of atomic properties. The results of such a calcu- 
lation for Ba are given in Table IIIII The one-particle e 
energies for ns and np orbitals in this calculation are set 



to ins = —0.48 a.u. and Cnpi/j 3/2 ~ —0.40 a.u., the ener- 
gies for the other partial waves are set to the Dirac-Fock 
values for the lowest orbital just as in the previous cal- 
culation. The results of this calculation are in excellent 
agreement with experiment. Moreover, the A('^£'i— ^Pq) 
energy difference which is very difficult to accurately cal- 
culate agrees with experiment to 11 cm^^. 

Our calculations have yielded accurate wave functions 



9 



TABLE III: Comparison of the CI + all-order results for the 
energy levels of Ba with experiment. Two-electron binding 
energies are given in the first row, the other values are counted 
from the ground state. The one-particle e energies for ns 
and np orbitals are set to e„s — —0.48 a.u. and 3/2 ~ 

—0.40 a.u., the energies for the other partial waves are set 
to the Dirac-Fock values for the lowest orbital. The energies 
are given in cm~^. The relative difference with experimental 
values is given in the last column in %. 



State 


Expt. 


Cl-f all-order 


A(%) 


6s' 'So 


122721 


122757 


-0.03 


6s5d^Di 


9034 


9012 


0.25 


6s5d^L>2 


9216 


9202 


0.14 




9597 


9603 


-0.07 


6s5d^D2 


11395 


11407 


-0.11 


Qs6p^Po 


12266 


12235 


0.26 


6s6p^Pi 


12637 


12602 


0.27 


6s6p^P2 


13515 


13491 


0.18 


6s6p^Pi 


18060 


18120 


-0.33 


5d6p^F2 


22065 


22201 


-0.62 


MGp^Fs 


22947 


23174 


-0.99 


MQp^Fi 


23757 


24005 


-1.04 


A(^Di - =*Po) 


3232 


3223 


0.28 



for subsequent evaluation of the atomic properties. Ma- 
trix elements of one-body operators such as El, E2, hy- 
perfine, parity-violation, etc., are evaluated in the frame- 
work of the CI -|- MBPT approach in the RPA approx- 
imation [1^ (as described in Step 6 of the previous sec- 
tion), sometimes with subsequent addition of the dom- 
inant normalization and structure radiation terms. In 
the CI + all-order method, we can use exactly the same 
method to evaluate matrix elements. Our preliminary 
calculations of the ^Pq polarizability values for Ca and 
Sr indicate better agreement of the CH-all-order ah ini- 
tio results with recommended values from Ref. [s^ in 
comparison with the CI-fMBPT approach. The com- 
plete implementation of the all-order approach would re- 
quire addition of the all-order corrections to the matrix 



elements beyond the modification of the wave function 
that is the subject of the present paper. This approach 
will implicitly include dominant normalization, structure 
radiation, and other corrections to all orders. In further 
work, we plan to replace the RPA matrix elements by all- 
order counterparts that are linear or quadratic functions 
of the excitation coefRcients 2] . The terms that arc ac- 
counted by the CI will have to be removed to avoid double 
counting. The ability to conduct calculations in various 
approximations will also allow one to carry out the eval- 
uation of uncertainties of atomic properties needed for 
many application, such as calculation of BBR shifts and 
the study of fundamental symmetries. 

VI. CONCLUSION 

We have developed a theoretical method combining 
the all-order approach currently used in precision cal- 
culations of properties of monovalent atoms with the 
configuration-interaction approach that is applicable for 
many-electron systems. This approach has been tested 
on the calculation of energy levels of divalent systems 
from Mg to Hg. We have demonstrated an improvement 
of at least a factor of three in agreement with experi- 
mental values for the two-electron binding energies and 
most excited state energies in comparison with the CI -I- 
MBPT method. Further work on this method will in- 
clude addition of the all-order terms beyond the RPA 
in the treatment of the transition amplitudes and other 
matrix elements for precision calculation of atomic prop- 
erties of systems with few valence electrons. 
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